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ABSTRACT 

We develop and apply an analytic method to predict peak counts in weak-lensing surveys. It is based on the theory of Gaussian 
random fields and suitable to quantify the level of spurious detections caused by chance projections of large-scale structures as well 
as the shape and shot noise contributed by the background galaxies. We compare our method to peak counts obtained from numerical 
ray-tracing simulations and find good agreement at the expected level. The number of peak detections depends substantially on the 
shape and size of the filter applied to the gravitational shear field. Our main results are that weak-lensing peak counts are dominated 
by spurious detections up to signal-to-noise ratios of 3-5 and that most filters yield only a few detections per square degree above this 
level, while a filter optimised for suppressing large-scale structure noise returns up to an order of magnitude more. 
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1. Introduction 

Wide-area surveys for weak gravitational lensing can be and 
have been used for counting peaks in the shear signal, which are 
commonly interpreted as the signature of sufficiently massive 
dark-matter halos. However, such detections are clearly contam- 
inated by spurious detections caused by the chance superposition 
of large-scale structures, and even by the shape- and shot-noise 
contributions from the background galaxies used to sample the 
foreground shear field. As a function of the peak height, what is 
the conttibution of genuine halos to these detections, and how 
much do the large-scale structure and the other sources of noise 
contribute? 

Given the pow er of lensing-peak number counts as a cos- 
mological probe jMarian et all 120091 iKratochvil et all 120091: 
iDietrich & Har tlap 2009), we address this question here after de- 
veloping a suitable analytic approach based on peak counts in 
Gaussian random fields, as laid out bv lBardeen etaT] dl986h . It 
is reasonable to do so even though at least the high peaks are 
caused by halos in the non-Gaussian tail of the density fluc- 
tuations because the noise and large-scale sttuctures contribu- 
tions remain Gaussian, and thus at least the contamination of the 
counts can be well described analytically. Peaks with the high- 
est signal-to-noise ratios are expected to be more abundant than 
predicted based on Gaussian random fields. 

Weak-lensing data are filtered to derive peak counts from 
them. Several linear filters have been proposed and used in the 
literature. They can all be seen as convolutions of the measured 
shear field with filter functions of different shape s. Many shapes 
have been proposed for different purposes ([Schneide r et alj 
U998tlSchirmer et alj|2004 iMaturi et alj|2005l) . One filter func- 
tion, called the optimal filter later on, was designed specifically 
to suppress the contribution from the large-scale structure by 
maximising the signal-to-noise ratio of halo detections against 
the shear field of the large-scale structure. 



We study three such filters here, with the optimal filter 
among them. Results will differ substantially, arguing for a care- 
ful filter choice if halo detections are the main goal of the ap- 
plication. We compare our analytic results to a numerical simu- 
lation and show that both agree at the expected level. We begin 
in § 2 with a brief summary of gravitational lensing as needed 
here and describe filtering methods in § 3. We present our an- 
alytic method in § 4 and compare it to numerical simulations 
in § 5, where we also show our main results. Conclusions are 
summarised in § 6. In Appendix A, we show predictions of peak 
counts and the noise levels in them for several planned and on- 
going weak-lensing surveys. 

2. Gravitational lensing 

Isolated lenses are characterised by their lensing potential 

W^^Trk f ®(D A 6,z)dz, (1) 
c z D d D s J 

where O is the Newtonian gravitational potential and Z\d,ds are 
the angular-diameter distances between the observer and the 
source, the observer and the lens, and the lens and the source, 
respectively. The potential tfr relates the angular positions /? of 
the source and of its image on the observer's sky through the 
lens equation 

p = e-Vxj,. (2) 

Since sources such as distant background galaxies are much 
smaller than the typical scale on which the lens properties 
change and the angles involved are small, it is possible to lin- 
earise Eq. © such that the induced image distortion is expressed 
by the Jacobian 

a .a / 1 - gi 82 \ 
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where k = V 2 ifr/2 is the convergence responsible for the isotropic 
magnification of an image relative to its source, and g(0) = 
y(0)/[l - k(0)] is the reduced shear quantifying its distortion. 
Here, 71 = (^,11 — ft, 22) /2 and 72 = <A,i2 are the two compo- 
nents of the complex shear. Since the angular size of the source GtanhW = (l 
is unknown, only the reduced shear can be estimated starting 
from the observed ellipticity of the background sources, 



(2) A filter optimised for halos with NFW density pro- 
file, approximating t heir shear signal with a hyperbolic tangent 
(ISchirmer et alj|2004 . 
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\ tanh (x/x c ) 



(9) 



e = 



e s +g 
1 + g*e s 



(4) 



where e s is the intrinsic ellipticity of the source and the asterisk 
denotes complex conjugation. 

3. Measuring weak gravitational lensing 

3.1. Weak lensing estimator 

In absence of intrinsic alignments bet ween background galax- 
ies due to possible tidal interactions dHeavens & Peac ock 1988; 
[Schneider & Bridle 2009), the intrinsic source ellipticities in 
Eq. dU average to zero in a sufficiently large source sample. An 
appropriate and convenient measure for the lensing signal is the 
weighted average over the tangential component of the shear y t 
relative to the position on the sky, 



AW 
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d 2 , 7t (0',e)Q(\e'-0\) . 



(5) 



The filter function Q determines the statistical properties of the 
estimator A. We shall consider three filter functions here which 
will be described in § 13.21 

Data on gravitational lensing by a mass concentration can be 
modeled by a signal s{0) = At{0) described by its amplitude A 
and its radial profile t, and a noise component n(0) with zero 
mean, e.g. 

y t (0) = At{0) + n(0) (6) 
for the tangential shear. The variance of the estimator A in (0 is 



(7) 



where W(k) is the frequency response of the survey depending on 
its geometrical properties, Q(k) is the Fourier transform of the 
filter Q, and P(k) is the power spectrum of the noise component. 

3.2. Weak lensing filters 

Different filter profiles have been proposed in the literature de- 
pending on their specific application in weak-lensing. We adopt 
three of them here which have been used so far to identify halo 
candidates through weak lensing. 

(1^ ) The polynomial filter described by ISchneider et alj 
(1199? 



6x 



TT V \ / 



(8) 



where the projected angular distance from the filter centre, x = 
r/r s , is expressed in units of the filter scale radius, r s , and H 
is the Heaviside step function. This filter was originally pro- 
posed for cosmic-shear analysis but several a uthors have used it 
also for dark matter h alo searches (see for e.g. lErben et alj2000t 
ISchirmer et aD2004l) . 



where the two exponentials in parentheses are cut-offs imposed 
at small and large radii (a = 6, b = 150, c = 50, and 
d = 47) and x c is a parameter defining the filter-profile slope. 
A good choice for the latt er is x c = 0.1 as empirically shown by 
iHetterscheidt et all (120051) . 

(3) The optimal linear filter introduced by iMaturi et al.l 
(2005) which, together with the optimisation with respect to the 
expected halo-lensing signal, optimally suppresses the contami- 
nation due to the line-of-sight projection of large-scale structures 
(LSS), 



Q opt (k) = a 



r(k) 
P(k) 



with 



P(k) 



(10) 



Here, f(k) is the Fourier transform of the expected shear profile 
of the halo and P(k) = P & + Plss(^) is the complete noise power 
spectrum including the LSS through Plss and the noise contri- 
butions from the intrinsic source ellipticities and the shot noise 
by P g = o- 2 g /(2n g ), given their angular number density n g and 
the intrinsic ellipticity dispersion <x g . This filter depends on pa- 
rameters determined by physical quantities such as the halo mass 
and redshift, the galaxy number density and the intrinsic elliptic- 
ity dispersion and not on an arbitrarily chosen scale which has 
to be dete rmined empirically throug h costly numerical simula- 
tions (e.g. lHennawi & S pergel 20 05b. An applica tion of this fil- 
ter to_flie_GaBoDS_survey (Schirmeret al. 2003) was presented 
in lMaturi et al.l d2007l), while a detailed comparison of these fil- 
ters was performed bv lPace et al. (2003) by means of numerical 
ray-tracing simulations. They found that the optimal linear filter 
given by Eq. ( fTOb returns the halo sample with the largest com- 
pleteness (100% for masses M > 3 x lO 14 ^ 1 M and ~ 50% for 
masses M ~ 2 x 10 14 h~ l M Q for sources at z s = 1) and the small- 
est number of spurious detections caused by the LSS (< 10% for 
a signal-to-noise threshold of S/N ~ 5). 

3.3. Weak lensing estimator and convergence 

In order to simplify comparisons with numerical simulations, we 
convert the estimator A from Eq. © to an estimator of the con- 
vergence, 



A(0) = j ^&K{ff)U{\ff -0|) , 
where U is related to Q by 

2 r e 

Q{6) = ^ defer u{&) - u(e) 



(ii) 



(12) 



dSchneiderlll996l) if the weight function U(ff) is defined to be 
compensated, i.e. 



AffffU(ff) = . 



(13) 



Equation dT2l has the form of a Volterra integral equation of 
the first kind which can be solved with respect to U once Q is 
specified. If lim x ^o Q(x)/x is finite, the solution is 



U(0) = -Q(6) ■ 



r ,2 , 

Jo v 



(14) 



Maturi et al.: Weak-lensing peak detections 



3 




x [arcmin] x [arcmin] 



Fig. 1. Overview of different weak-lensing filters. The left panel shows the three filters adopted here to be used on shear catalogues, 
while the central and right panels show the corresponding filters to be used on convergence fields both in real and Fourier space, 
respectively. For illustration only, the spatial frequencies are rescaled such that the main filters peaks coincide. 



dPolyanin & Manzhirovl [l998) which can be solved analytically 
for the polynomial filter 

t/polyW = ^ (l - X 2 )(^ - * 2 )//(l - X) , (15) 

and numerically for the hyperbolic-tangent filter of Eq. © 
with an efficient recursive scheme over the desired radii 9. If 
lim x ^o Q(x)/x = °o as in the case of the optimal filter, Eq. (fT4] i 
can be solved by introducing an exponential cut-off at small radii 
to avoid the divergence. The correct solution is obtained if the 
cut-off scale is close to the mean separation between the back- 
ground galaxies, so that no information is lost. Alternatively, 
Eq. ( fTSl i can be solved iteratively with respect to Q by 



U (S) = -Q(0) , U„(0) = -Q(6) + 



2 f 6 
~& Jo 



AffffU n ^{&) . (16) 



The iterative procedure is stopped once the difference U„(0) - 
U„-\{6) is sufficiently small. After U(0) has been found, an ap- 
propriate constant c has to be added in order to satisfy the com- 
pensation requirement, Eq. ( 113) . It is given by 



"mas 

— f 

p max +J 



dffffU{ff) 



(17) 



We show in Fig. (Q3 the resulting filter profiles to be used on 
shear catalogues through Eq. (O and their corresponding vari- 
ants to be used on convergence fields with Eq. (fTTT i both in 
real and Fourier space. All of them are band-pass filters and 
the two of them designed for halo searches have larger ampli- 
tu des at higher frequenc ies compared to the polynomial filter 
bv lSchneider et al.l (119981) . where the halo signal is most signifi- 
cant. This feature is particularly prominent for the optimal filter, 
which is additionally negative at low frequencies, where the LSS 
signal dominates. These two features ensure the minimisation of 
the LSS contamination in halo searches. 



4. Predicting weak lensing peak counts 

Our analytic predictions for the number counts of weak-lensing 
detections as a function of their signal-to-noise ratio are based on 
modelling the analysed lensing data, resulting from Eq. ( fTTl i. as 
an isotropic and homogeneous Gaussian random field. This is an 
extremely good approximation for the noise and the LSS com- 
ponents, but not for the non-linear structures such as sufficiently 
massive halos, as we shall discuss in Sec. 



4.1. Statistics of Gaussian random fields 

An n- dimensional random field F(r) assigns a set of random 
numbers to each point r in an n-dimensional space. A joint prob- 
ability function can be declared for m arbitrary points rj as the 
probability to have field values between F{rj) and F(rj)+dF(r j), 
with j — 1 , . . . , m. For Gaussian random fields, the field it- 
self, its derivatives, integrals and any linear combination thereof 
are Gaussian random variables which we denote by y, with 
mean value (y,) and central deviations Ay, := y, ■ - (y,), with 
i = I,..., p. Their joint probability function is a multivariate 
Gaussian, 



f > (y u ...,y v )dyi---dy p = 



with the quadratic form 



1 



e^dyi — dyfe (18) 



yj{2nf det(Af) 



fl:=l|>(Al %A,. 



(19) 



where M is their covariance matrix with elements Ai/j : = 
(Ay/Ayj). All statistical properties of a homogeneous Gaussian 
random field with zero mean are fully characterised by the two- 
point correlation function £(ri,r2) = £(\r\ — r%\) '■= (F(r{)F(r2)) 
or equivalently its Fourier transform, the power spectrum P(k). 
In our case, this is the sum of the matter fluctuation convergence 
power spectrum, PlssW, an d observational noise, P n (k). 

Since we are interested in gravitational-lensing quantities 
such as the convergence k, we here consider two-dimensional 
Gaussian ran d om fi elds only. We adopt the formalism of 
Bard een et al.l d 1986b . where F — k, 77, = djF and £y = djdjF 
denote the convergence field and its first and second derivatives, 
respectively. 



4.2. Definition of detections: a new up-crossing criterion 

We define as detection any contiguous area of the field k which 
exceeds a given threshold, — S/N • cr^ , determined by the re- 
quired signal-to-noise ratio, S/N, and the variance cr^ of an esti- 
mator A (see Eq. [7J. This definition is widely used in surveys for 
galaxy clusters or peak counts in weak-lensing surveys and can 
easily be applied both to real data and Gaussian random fields. 

Each detection is delimited by its contour at the threshold 
level Kfa. If this contour is approximately circular, it has a single 
point r up , called up-crossing point, where the field is rising along 
the x-axis direction only, i.e. where the field gradient has one 
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vanishing and one positive component (see the sketch for type-0 
detections in the lower panel of Fig. |2j, 



F(r ap ) = Kth , m( r w) = o . m( r u P ) > o 



(20) 



Since we assume k to be a homogeneous and isotropic random 
field, the orientation of the coordinate frame is arbitrary and ir- 
relevant. The conditions expressed by Eq. < f2Qb define the so- 
called up-crossing criterion which allows to identify the detec- 
tions and to derive their statistical properties, such as their num- 
ber counts, by associating their definition to the Gaussian ran- 
dom field variables F, tj\ and 772. 

However, this criterion is prone to fail for low thresholds, 
where detections tend to merge and the isocontours tend to devi- 
ate from the assumed approximately circular shape. This causes 
detection numbers to be overestimated at low cut-offs because 
each "peninsula" and "bay" of their profile (see type-1 in Fig. [2j 
would be counted as one detection. We solve this problem by 
dividing the up-crossing points into those with positive (red cir- 
cles) and those with negative (blue squares) curvature, f n > 
and < respectively. In fact, for each detection, their dif- 
ference is one (type-1) providing the correct number count. The 
only exception is for those detections containing one or more 
"la goons" (type-2) since each of them decreases the detection 
count by one. But since this is not a frequent case and occurs 
only at very low cut-off levels, we do not consider this case here. 

4.3. The number density of detections 

Once the relation between the detections and the Gaussian ran- 
dom variables y = {k±, r\\, 772, fn) and their constraints from 
Eq. ( |20l together with > or (n < are defined, we 
can describe their statistical properties through the multivariate 
Gaussian probability distribution given by Eq. ( fTST l with the co- 
variance matrix 
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3o-2/8 



(21) 



Thi s matrix differs from that deri ved by iBardeen et alJ (l986) 
and lAngrick & Bartelmannl (12009) because we are dealing with 
a 2-dimensional rather than a 3-dimensional field. Here, the <r ; 
are the moments of the power spectrum P(k) 



Fig. 2. Weak lensing detection maps. The top four panels show 
the segmentation of a realistic weak-lensing S/N map for in- 
creasing thresholds: 0.1, 0.5, 1, and 2, respectively. The bot- 
tom panel sketches the three discussed detections types together 
with the points identified by the standard and the modified up- 
crossing criteria. Red circles and blue squares correspond to up- 
crossing points for which the second field derivatives are (\ \ > 
and £n < 0, respectively. 



•5-J 



k 2j+1 dk^ 
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-P(k)W 2 (k)\Q(k)\ 2 



(22) 



so that the infinitesimal volume element d/cd^i can be written as 
dxdrji = I det J|d 2 r, where J is the Jacobian matrix, 



where W(k) is the frequency response of the survey given by 
its geometry (see Sec. I5.21 > and Q(k) is the Fourier transform of 
the filter adopted for the weak lensing estimator (see Sec. 13.21 ). 
The determinant of M is (3o- 2 fr A l o-\ - 2o"j)/32 and Eq. ( fT9l ) can 
explicitly be written as 



Q 



'2if Hti°o + 8fiiKo-f +3kV; 



3o >2 



2a\ 



(23) 



Both k and rj\ can be expanded into Taylor series around the 
points r up where the up-crossing conditions are fulfilled, 



k (r) ~ Kth + 2] m( r - . m Or) ~ ^ i u(* 



'"up)/ 



(24) 
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drill dx\ dr]\l dx 2 



m m 
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(25) 



and I det J I = 1 772^1 1 1 since 771 = 0. The number density of up- 
crossing points at the threshold k± with > 0, and < 0, n + 
and n~ respectively, can thus be evaluated as 



^(/Cth) = ± 



j d m J d£ n \m£n\'P(K = i«h,m = 0,ift,?n) ,(26) 




where V(k, 771,772,^11) is the multivariate Gaussian defined by 
Eq. dT~8b with p - 4, the correlation matrix ( |2TI ). and the 
quadratic form d23l . Both expressions can be integrated analyt- 
ically and their difference, ndet^th) = « + (Kth) - « (*th), as ex- 
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plained in Sect. 14.21 returns the number density of detections 
«det above the threshold 



«<kt(*th) = — p— ~r — — exp 
4 V2tt 3 / 2 Wo o- () 



20%, 



(27) 



Note how the dependence on cr 2 drops out of the difference n + - 
n~, leading to a very simple result. 

For completeness we report the number density estimate also 
for the classical up-crossing criterion, Eq. ( l20t only, where the 
constraint on the second derivative of the field, £\ \ , is not used, 



«upOth) 



1 



4 V2tp 



o-o 



— exp 



2o% 



(28) 
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with y :- ^cr^cr^ - 2cr|. This number density converges to 
the correct value njet for k± — > oo, i.e. large thresholds, because 
erf(jc) — > 1 and exp(-x 2 ) — > 0. This reflects the fact that, for 
large thresholds, the detection shapes become fully convex and 
any issues with more complex shapes disappear. 



5. Analytic predictions vs. numerical simulations 

We now compare the number counts of detections predicted by 
our analytic approach with those resulting form the analysis of 
synthetic galaxy catalogues produced with numerical ray-tracing 
simulations. 



5.1. Numerical simulations 

We use a hydrodynamical, nu merical jV-body simulation carried 
out with the code GADGET-2 dSpringelll2005l) . W e briefly sum- 
marise its main characteristics here and refer to iBorgani et al.l 
(I2004I) for a more detailed discussion. The simulation represents 
a concordance ACDM model, with dark-energy, dark-matter and 
baryon density parameters iX\ = 0.7, Q m = 0.3 and £2b = 0.04, 
respectively. The Hubble constant is Hq = lOO/ikms -1 Mpc~' 
with h = 0.7, and the linear power spectrum of the matter-density 
fluctuations is normalised to cr 8 = 0.8. The simulated box is a 
cube with a side length of 192 /i^'Mpc, containing 480 3 dark- 
matter particles with a mass of 6.6 x 10 9 h~ l M Q each and an 
equal number of gas particles with 8.9 x 10 8 h~ l M & each. Thus, 
halos of mass 10 13 hr x M & are resolved into several thousands of 
particles. The physics of the gas component includes radiative 
cooling, star formation and supernova feedback, assuming zero 
metallicity. 

This simulation is used to construct backward light cones by 
stacking the output snapshots from z = 1 to z = 0. Since the 
snapshots contain the same cosmic structures at different evo- 
lutionary stages, they are randomly shifted and rotated to avoid 
repetitions of the same cosmic structures along one line-of-sight. 
The light cone is then sliced into thick planes, whose parti- 
cles are subse quently projected with a tria ngular-shaped-cloud 
scheme (TSC, iHocknev & Eastwoodl H988) on lens planes per- 
pendicular to the line-of-sight. We trace a bundle of 2048 x 2048 
light rays through the light cone which start propagating at the 
observer into directions on a regular grid of 4.9 degrees on each 
side. The effective resol ution of this ray -tracing simulation is of 
1' (for further detail, see Pace et al. 20071). 

The effective convergence and shear maps obtained from the 
ray-tracing simulations are used to lens a background source 



le-08 
1e-09 
1e-10 
1 e-1 1 
1e-12 
1e-13 
1e-14 





: V ] I 




simulation 


'II. : 


■ *(P L SS+Pn) — — - , 





100 



1000 



10000 




100 



1000 



10000 



1e-07 



1e-08 



1e-09 



1e-10 



1 e-1 1 



1e-12 




simulation 
*(PLSS+Pn) 



100 



1000 



10000 



Fig. 3. Comparison between the power spectrum measured from 
the analysis of shear catalogues derived from numerical sim- 
ulations (heavy line) a nd the prediction resulting from the 
Peac ock & Doddsl (1 19961) power spectrum convolved with the 
survey frequency response and the adopted filter. We show from 
top to bottom the results for the polynomial filter with r s = 5'5, 
the hyperbolic-tangent filter with r s = 10' and the optimal filter 
with the cluster model scale radius set to r s - 1'. 

population according to Eq. (0). Galaxies are randomly dis- 
tributed on the lens plane at z — 1 with a number density of 
n„ - 30 arcmin -2 and have intrinsic random ellipticities drawn 
from the distribution 



Pied = 



exp|(l-ksl 2 )A^| 
7rcr 2 |exp(l/cT2)-l| 

where o> = 0.25. 



(29) 
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Fig. 4. Number of peaks detected in the numerical simulation 
(shaded area) compared to the prediction obtained with the 
proposed method both with the original up-crossing criterion 
(dashed line) and with the new blended up-crossing criterion 
(points with error bars). The standard up-crossing criterion is 
a good approximation for high signal-to-noise ratios but fails for 
lower S/N, which are well described by the new version. 



Synthetic galaxy catalogues produced in this way are finally 
analysed with the aperture mass (Eq. [5]) evaluated on a regular 
grid of512x512 positions covering the entire field-of-view of 
the light cone. All three filters presented in § |3.2| were used with 
three different scales: the polynomial filter with r s - 2.'75, 5'.5, 
and 11', the hyperbolic-tangent filter with r s = 5', 10', and 20', 
and the optimal filter with scale radii of the cluster model set to 
r s = 1', 2', and 4'. These scales are chosen such that the effective 
scales of the different filters can be compared. For a statistical 
analysis of the weak-lensing detecti ons and their rela tion to the 
numerical simulations structures, see iPace et all d2007l) . 



5.2. Accounting for the geometry of surveys: the window 
function 

Our analytic predictions for the number density of detections 
accounts for the full survey geometry, through the frequency re- 
sponse W(k) in Eq. ( 1221 . Any survey geometry can be consid- 
ered but, for sake of simplicity, we evaluate our estimates for a 
square-shaped field-of-view without gaps. Under these assump- 
tions, the frequency response, W(k), is the product of a high- 
pass filter suppressing the scales larger than the light cone's side 
length Lf = 2n/kf = 4.9 deg, 



W{(k) = exp 



SI 

k 2 ) 



(30) 



a low-pass filter imposed by the average separation d = 27i/k s = 



-1/2 

1 ' 



0'18 between the galaxies, 



W 2 Jk) = exp 



k 2 



(31) 



and a low-pass filter related to the resolution R plx = 0'.57 used to 
sample the sky with the weak lensing estimator A of Eq. (0, 



2 yfn IkR 



kR 



pix 



v pix 



(32) 



where i\(x) is the cylindrical Bessel function of order one. The 
latter function is a circular step function covering the same area 
as a square-shaped pixel of size R p \ x . The square shapes of the 
field-of-view and the pixels could be better represented by the 
product of two step functions in both the x- and y-direction, but 
the low gain in accuracy does not justify the higher computa- 
tional cost. Finally, for the comparison with our numerical ray- 
tracing simulation, we have to account for its resolution proper- 
ties which act on the convergence power spectrum only by con- 
volving P K with a low-pass filter 



W 2 (k) 



' k. 



(33) 



where k s = 2n/l arcmin 1 as discussed in § 15.11 

Note that when relating the detection threshold to the signal- 
to-noise ratio S/N according to the estimator variance using 
Eq. (0 and = S/N • crj, all filters mentioned are used except 
for Wpi x , which, of course, does not affect the estimator variance. 



5.3. Comparison with numerical simulations 

We first compare in Fig.|3]the expected filtered power spectra for 
the convergence and the intrinsic noise of the galaxies, P|WQ| 2 , 
assumed in Eq. (1221 1. with those measured from the analysed 
synthetic galaxy catalogue. All main features are very well re- 
produced. Only at high frequencies the assumed power spectra 
drop slightly more steeply than measured in the numerical sim- 
ulations. 

A comparison of the original up-crossing criterion with the 
new blended up-crossing criterion presented here is shown in 
Fig. [4] together with the number counts obtained from the nu- 
merical simulations. Only the result for the optimal filter with 
r s = 1' is shown for clarity. As expected, the two criteria agree 
very well for high signal-to-noise ratios since the detections are 
mostly of type-0, i.e. approximately circular, as shown in the left 
panel of Fig. [2] while the merging of detections at lower signal- 
to-noise ratios is correctly taken into account only by our new 
criterion. Note that our analytic approach approximates the data 
as random Gaussian fields, very well representing noise and LSS 
but not non-linear objects such as galaxy clusters. Thus, partic- 
ular care has to be taken when comparing the predicted num- 
ber counts with real or simulated data by modelling the non- 
linear structures, which is difficult and uncertain, or by avoiding 
their contribution in first place. We follow the latter approach by 
counting the negative instead of the positive peaks found in the 
galaxy catalogues. In fact, massive halos contribute only posi- 
tive detections in contrast to the LSS and other sources of noise 
which equally produce positive and negative detections with the 
same statistical properties. 

This is an interesting application for cosmology, where a ro- 
bust prediction of number counts in the linear regime only can 
be directly compared to data or, turning the argument around, 
can be used to statistically correct the halo number counts by 
using the data only. In fact, the difference of the positive and 
negative detection counts is caused by non linear structures only 
and their Poisson fluctuations. 

Our analytic predictions for all filters and both positive and 
negative detection counts resulting from the synthetic galaxies 
catalogue from the numerical simulation are shown in Fig. [5] 
The high signal-to-noise ratio tail caused by the nonlinear struc- 
tures is present only in the positive detection counts, as expected. 
The agreement with the negative detections is within the 1-cr 
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Fig. 5. Nu mber of weak tensing pe aks, shown as a func t ion of the sign al-to-noise ra t io, p redicted with the analytic method presented 
here for the I Schneider et all d!998l) . the lSchirmer et alj d2004l) and the M aturi et"al] d2005l) filters from top to bottom, and increasing 
filter radii from left to right. The number counts generated by the intrinsic galaxy noise alone, P g , and the LSS alone, Plss, are also 
shown. Numbers refer to a survey of one square degree with a galaxy number density of n g = 30 arcmin -1 and an intrinsic shear 
dispersion of cr g = 0.25. The results are compared with the number counts of positive as well as negative peaks detected based on 
the synthetic galaxy catalogues from the numerical simulation. 



error bars (considering the uncertainties for a one square de- 
gree survey) except for lSchirmer et alj d2004l) filter and the op- 
timal filter, both with a scale of 5 arcmin, which are compatible 
only at a 2-<x level. It is plausible that these small deviations are 
caused by the ignored correlation between the pixels which, in 
the noise evaluation, are treated as independent. The filter radii 
are in fact larger than the pixel separation, so that different pixels 
have some common galaxies. 

We finally compare the contribution of the LSS and the noise 
to the total signal by treating them separately. Their number 
counts are plotted with dashed lines in Fig. [5] All filters show an 
unsurprisingly large number of detections caused by the noise up 
to signal-to-noise rations of 3 and a number of detections caused 
by the LSS increasing with the filter scale except for the optimal 
filter, which always suppresses their contribution to a negligible 
level. Thus, the LSS contaminates halo catalogues selected by 
weak lensing up to signal-to-noise ratios of 4 - 5 if its contribu- 
tion is ignored in the filter definition. Note that the total number 
of detections can be obtained only by counting the peaks from 
the total signal, i.e. LSS plus noise, and not by adding the peaks 
found in the two components separately, because the blending of 
peaks is different for the two cases. 



6. Conclusion 

We have developed an analytic method for predicting peak 
counts in weak -lensing surveys, ba sed on the theory of Gaussian 
random fields (Bar deen et al.f l986). Peaks are typically detected 
in shear fields after convolving them with filters of different 
shapes and widths. We have taken these into account by first 
filtering the assumed Gaussian random field appropriately and 
then searching for suitably defined peaks. On the way, we have 
argued for a refinement of the up-crossing criterion for peak de- 
tection which avoids biased counts of detections with low signal- 
to-noise ratio, and implemented it in the analytic peak-count pre- 
diction. Peaks in the non-linear tail of the shear distribution are 
underrepresented in this approach because they are highly non- 
Gaussian, but our method is well applicable to the prediction of 
spurious counts, and therefore to the quantification of the back- 
ground in attempts to measure number densities of dark-matter 
halos. We have compared our analytic prediction to peak counts 
in numerically simulated, synthetic shear catalogues and found 
agreement at the expected level. 

Our main results can be summarised as follows: 

- The shape and size of the filter applied to the shear field has a 
large influence on the contamination by spurious detections. 
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For the optimal filter, the contribution by large-scale struc- 
tures is low on all filter scales, while they typically contribute 
substantially for other filters. 

- Shape and shot noise due to the background galaxies used to 
measure the shear from are a large source of spurious peak 
counts for all filters, and the dominant source for the optimal 
filter. 

- Taken together, large-scale structure and galaxy noise con- 
tribute the majority of detections up to signal-to-noise ratios 
between 3-5. Only above this level, detections due to real 
dark-matter halos begin dominating. 

- The optimal filter allows the detection of ~ 30-40 halos per 
square degree at signal-to-noise ratios high enough for sup- 
pressing all noise contributions. For the other filters, this 
number is lower by almost an order of magnitude. 

Our conclusions are thus surprisingly drastic: peak counts in 
weak-lensing surveys are almost exclusively caused by chance 
projections in the large-scale structure and by galaxy shape and 
shot noise unless only peaks with high signal-to-noise ratios are 
counted. With typical filters, only a few detections per square de- 
gree can be expected at that level, while the optimal filter returns 
up to an order of magnitude more. 
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which cannot be considered in this paper in any detail. Finally we fixed the sam- 
pling scale, described by Eq. (32), to be 5 times smaller than the typical filter 
scale in order to avoid undersampling, i.e. such that the high frequency cut-off is 
imposed by the filters themselves. The results are shown in Tab. [Altogether with 
the number counts obtained with a simple Gaussian fi lter, usually used to gether 
with the Kaiser & Squires shear inversion algorithm (Kaiser & Squires 1993). 
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Appendix A: Forecast for different weak lensing 
surveys 

For convenience, we evaluate here the expected number density of peak counts 
for a collection of present and future weak-lensing surveys. To give typical val- 
ues, we assumed for all of them a square-shaped field of view, a uniform galaxy 
number density and no gaps for two main reasons. First, their fields-of-view are 
typically very large and thus do not affect the frequencies relevant for our evalu- 
ation. Second, the masking of bright objects can be done in many different ways 
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Table A.l. Expected number counts of peak detections per square degree for S/N = 1,3,5 and for several present and future 
weak-lensing surveys with different intrinsic ellipticity dispersion, cr e , and galaxy number density, n g , per arcmin 2 . For each filter, 
we used three different scales, namely Q po i y : scale- 1=2.'75, scale-2=5.'5, scale-3=ll'; Q tan h-' scale-l=5', scale-2=10', scale-3=20'; 
Q opt : scale-l = 10 14 M o and scale-2=5 x 1O 14 M . Q gaus (Gaussian FWHM): scale-l=l', scale-2=2', scale-3=5'. 



Pan-STARRS 

cr e = 0.3, n g = 5 


Qpoly 

1 3 5 


Qtanh 

1 3 5 


Qopt 

1 3 5 


Qgaus 

1 3 5 


scale- 1 


207.7 8.127 0.002 


252.7 8.824 0.002 


186.3 6.243 0.001 


3125 131.1 0.042 


scale-2 


51.49 2.82 0.001 


61.8 3.214 0.002 


62.69 3.576 0.002 


989.6 38.93 0.011 


scale-3 


12.45 1.258 0.002 


14.02 1.518 0.003 




173.2 7.82 0.003 



DES 

a € = 0.3, w g = 10 


Qpoh 

1 3 5 


Qtanh 

1 3 5 


Qopt 

1 3 5 


Qgaus 

1 3 5 


scale- 1 


206.6 9.55 0.004 


248.9 10.98 0.004 


288.8 12.72 0.005 


3593 144.5 0.043 


scale-2 


50.09 4.178 0.005 


56.87 5.112 0.008 


95.6 8.325 0.012 


1047 41.67 0.012 


scale-3 


11.67 2.339 0.0174 


11.92 2.847 0.030 




169.8 9.807 0.006 



CFHTLS 

cr € = 0.3, n g = 20 


Qpoh- 

1 3 5 


Qtanh 

1 3 5 


Qopt 

1 3 5 


Qgaus 

1 3 5 


scale- 1 


206.6 9.907 0.004 


249.6 11.48 0.004 


324 14.12 0.005 


3971 151.6 0.041 


scale-2 


49.76 4.545 0.007 


55.86 5.61 0.010 


104.5 9.519 0.015 


1085 42.15 0.012 


scale-3 


11.49 2.622 0.025 


11.51 3.166 0.044 




169.7 10.28 0.007 



Subaru 

cr € = 0.3, n g = 30 


Qpoh 

1 3 5 


Qtanh 

1 3 5 


Qopt 

1 3 5 


Qgaus 

1 3 5 


scale- 1 


198.5 16.22 0.020 


219.5 21.51 0.038 


603.2 42.99 0.04045 


4110 160.9 0.046 


scale-2 


44.84 10.82 0.117 


42.97 13.32 0.237 


172.8 29.81 0.1642 


1070 50.42 0.021 


scale-3 


9.406 6.321 0.528 


7.857 6.457 0.807 




151.6 19.14 0.057 



EUCLID 

cr £ = 0.3, n g = 40 


Qpoh 

1 3 5 


Qtanh 

1 3 5 


Qopt 

1 3 5 


Qgaus 

1 3 5 


scale- 1 


194.3 20.01 0.039 


206.3 27.29 0.088 


730.9 58.61 0.070 


4189 165.5 0.048 


scale-2 


42.64 14.25 0.295 


38.54 16.8 0.591 


197.7 40.75 0.321 


1062 54.99 0.027 


scale-3 


8.653 7.642 1.104 


6.873 7.282 1.514 




143.6 24.19 0.127 



LSST 

cr e = 0.22, n g = 50 


Qpoh- 

1 3 5 


Qtanh 

1 3 5 


Qopt 

1 3 5 


Qgaus 

1 3 5 


scale- 1 


174.8 42.42 0.463 


156.8 56.13 1.333 


1206 138 0.334 


4169 187.7 0.070 


scale-2 


34.5 28.13 3.464 


26.15 26.81 5.218 


269.2 95.09 2.198 


991.5 82 0.104 


scale-3 


6.403 10.32 4.964 


4.519 8.139 4.889 




113.5 48.96 1.688 



SNAP 

a e = 0.3, n g = 100 


Qpoh- 

1 3 5 


Qtanh 

1 3 5 


Qopt 

1 3 5 


Qgaus 

1 3 5 


scale- 1 


172.6 45.42 0.5824 


152.5 59.33 1.664 


1322 148.6 0.3481 


4287 190.2 0.069 


scale-2 


33.73 29.39 4.133 


25.22 27.43 6.009 


281.6 102.2 2.494 


991.3 85.32 0.117 


scale-3 


6.218 10.41 5.403 


4.355 8.1 5.19 




110.8 51.68 2.083 



